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INTRODUCTION 



The following is a brief description of all the relocatable math library 
subroutines available on tape number 099-000001. These subroutines appear in 
alphabetical order according to function. Each description will elaborate on the 
items explained below. 

PURPOSE: Explains the function performed by the routine. 



TTTT.R- 



ENTRY: 



I, i t ro o t-T-» a r> o *-*-» *-\ *-» -P *-"U f~, -*+*■ 



»_o u «, i.am t ui me iuuLiiic ^necessary ior eaiting me library tape). 

Gives the name by which a routine is referenced in an . EXTN 
statement. This name is identical to the JSR entry point unless 
specified otherwise. 



INPUT: Describes necessary input format. 

OUTPUT: Describes results of the routine. 

CALLING SEQUENCE AND ENTRY POINTS: 

The relocatable routines contained in the math library 
are called by declaring the appropriate entry point as 
a normal external within the user program. For example, 
to call double precision absolute value: 



EXTN 



.DABS 



JSR 



§DUMMY 



DUMMY: 



.DABS 



The names of the entry point(s) are given for all routines, and unless 
otherwise noted all user calls must use the above method. 

ERROR CONDITIONS: 

Explains or cautions about the idiosyncrasies of the routine. 



CARRY AND REGISTERS: 

Gives states of active registers upon exit. 

LENGTH AND TIME: 

Gives the number of words used by this routine and an approxi- 
mation of execution time. Unless otherwise noted, execution 
times are calculated for the Nova. To obtain approximate exe- 
cution times for other Nova -family machines, use the following 
conversions: 

. 2 * Nova execution time •■- Supernova or Nova 800 execution times 
. 35 * Nova execution time -.-. Nova 1200 execution time. 

ALGORITHM: 

Describes the method used to produce the desired result. 

REFERENCE: 

Cites literature that may be of use in obtaining further information. 

PROGRAM LISTING: 

An assembly listing of the routine is given. 

In most instances there will be no ENTRY, ALGORITHM, or PROGRAM 
LISTING entries given. The ALGORITHM and PROGRAM LISTING 
entries will be included in future editions of this manual. ENTRY infor- 
mation is given only when the entry is different from the JSR entry point. 
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in 



ARQOT TTTP WAT rTTJ 

(Single Precision) 



PURPOSE: 



TITLE: 



INPUT: 



OUTPUT: 



This routine computes the absolute value of a fixed point, 
single precision, two's complement number. 



The title is . ABS. 



The input is a single precision number in ACO. 



The absolute value of the input is returned in ACO. 



CALLING SEQUENCE AND ENTRY POINT: 

Indirect to . ABS with normal return to the instruction 
following the call. 

ERROR CONDITIONS: 

The absolute value of -2**15 cannot be represented and 
will be returned unchanged. 

CARRY AND REGISTERS: 

Carry and ACO may be destroyed; AC1, AC2 and AC3 are 
unchanged. 

LENGTH AND TIME: 

This routine consists of 3 words and is normally relocatable. 
For X » 0, execution time is 8.2 us. 
For X -= 0, execution time is 13. 8 us. 
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ABSOLUTE VALUE 
(Double Precision) 



PURPOSE: 



TITLE: 



INPUT: 



OUTPUT: 



This routine computes the absolute value of a double precision, 
fixed point, two's complement number. 



The title is . DABS. 



A number in ACO (high order), AC1 (low order). 



The absolute value of the input returned in AC 0, AC1 - high 
order in ACO, low order in AC1. 

CALLING SEQUENCE AND ENTRY POINT: 

Indirect at . DABS, with normal return to the instruction 
following the call. 

ERROR CONDITIONS: 

Caution: The absolute value of -2**31 cannot be represented 
and is returned unchanged. 

CARRY AND REGISTERS: 

ACO, AC1, and Carry are destroyed; AC2 and AC3 remain 
unchanged. 

LENGTH AND TIME: 

This routine consists of 6 instructions and is normally 

relocatable. 

For X =■ 0, execution is 8.2. us 

For X «= 0, execution is 19.4 us 
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PURPOSE: 



TITLE: 



INPUT: 



OUTPUT: 



ADDITION 



(Double Precision) 



This routine computes the sum of two double precision, two's 
complement integers. 



The title is . DADD. 



The first operand must be in ACO, AC1 (high order, low order). 
The second operand must be in storage, higher order word 
followed by lower order word. The address nf t-hp hio-hPi- 
order word of the second operand must be given after the 
JSR ©DUMMY. 



The double nrerisinn sum will J->o rotnimo^ i-n Apn \ ^ /u;~u 
order, low order). 



CALLING SEQUENCE AND ENTRY POINT: 

Indirect at . DADD, then address of second operand with 
return to the instruction following the second operand 
address. 

ERROR CONDITIONS: 

Caution: No check is made for overflow. 

CARRY AND REGISTERS: 

ACO, AC1, AC3 and Carry are destroyed; AC2 remains 
unchanged. 

LENGTH AND TIME: 

This routine consists of 15 (octal) instructions and is 
normally relocatable. 
Execution time is 54. 9 us. 
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ARCTANGENT 



(Single Precision) 

PURPOSE: 

To calculate the fixed point arctangent of the quotient 
of two input arguments. 



TITLE: 



ENTRY: 



ATANX 



.ATANX 



INPUT: 



OUTPUT: 



Argument dividend in ACO, argument divisor in AC1. 
Both arguments are expressed in radians in the following 



sign| integer 



fraction 



bit bits 1 and 2 bits 3 through 15 

The sign bit is set to a 1 only if the argument is negative. 

The result x, expressed in radians, falls in the range 

-TT <. x <TT . The result, in AC2, is given in the same format 
as that described for input arguments. 



CALLING SEQUENCE: 

JSR indirect through page zero entry . ATANX . Return 
is to the next sequential location following the call. 

ERROR CONDITIONS: 

None; all input arguments will be interpreted in the 
format illustrated above. 

CARRY AND REGISTERS: 

ACQ and AC1 are saved; AC3 and Carry are destroyed. 

LENGTH AND TIME: 

One ZREL and 100 octal NREL locations. Average 
execution time on the NOVA 1200 is 1. 3 ms. 



ARCTANGENT (con* t) 



(Single Precision) 
ALGORITHM: 

The quotient of the input arguments, x, is found by means 
of a call to the unsigned integer divide routine, DVD. For 
the range « x « 1, the arctangent of x is calculated to be 
equal to x * P (x**2). Calls to the unsigned integers multiply 
routine (.MPYU) and Polynomial expansion function (. POLY) 
a^made. A sixth order polynomial is computed, 
P(x ) = P Q + Pj x + P 2 x . . . + P 6 x 6 ,with the following 
coefficients: 



Pn = 


. 99999 


p 1 = 


- .33326 


p 2 = 


.19881 


p 3 = 


- .13487 


P 4 = 

p. = 


.83871 * 10" 1 
,37012 * 10" l 


p 6 = 


.78633* 10" 2 



For other values of x, one of the following quadrant 
adjustments is made, where m and n represent the original 
arguments input in ACO and AC1 respectively. 



ARCTAN ( —) =77/2 - ARCTAN ( | ) 

-m m 

ARCTAN C^— ) = -ARCTAN ( — ) 

+m 

ARCTAN (_ n ) = IX- ARCTAN (f") 



-m 



m 



ARCTAN <— )= - (TT- ARCTAN (f-) ) 

REFERENCE: 

John F. Hart, 'Computer Approximations ,, New York: John 
Wiley &Sons, Inc., 1968; pages 128 - 129, INDEX 4990 
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ARCTANGENT (cont'd) 
(Single Precision) 



PROGRAM LISTING: 









y 


FIXED POINT 








.TITLE 


ATANX 








. EN 1 


.ATAMX 








. EXTD 


.POLY* MPYU 








. EXTN 


DVD 








. r. KFL 




00000- 


000007 * 


. ATANX: 


ATANX 
.NREL 




« 


000000 


RTURN : 







CT0001 • 


000 000 


S A vest 


O 




raping? ' 


000000 


SAV1 : 







00 (7! (71 3 * 


000000 


SAV2: 







00004 ' 


000000 


SIGN: 







00^0 5' 


000000 


COM PL: 







00006" 


000000 


SUP pl: 







00007 


•054771 


ATAMX: 


ST A 


3, RTURN 


00(7! 1 Pj 


•040771 




ST A 


O,SAV0 


000 1 1 


•044771 




ST A 


1,SAV1 


H 1 P 


• 101 1P0 




M0V2L 


0,0 


000 13 


' 1 5P5 60 




SUBCL 


P, P 


n i /i 


•050770 




STA 


P, SIGN 


1 5 


• 1 P51P0 




MO V?.L 


1* 1 


0001 6 


• 15P5 60 




SUBCL 


P, 2 


0001 7 


•050767 




STA 


P, SUP PL 


00"! P0 


' 17 6400 




SUB 


3,3 


000 PI 


» IP! 643 P 




SUBS* 


0, 1,S£C 


000 PP 


•000 404 




JMP 


.+4 


00P3 


' 1 1 1000 




MOV 


0,P 


000P4 


•175400 




IMC 


3,3 


000 P 5 


•1 PI 001 




MOV 


1,0, SKP 


000P6 


•131000 




MOV 


1,P 


000P7 


•054756 




STA 


3, COM PL 


00030 


•11?/' 15 




SUB* 


0,P, SMH 


00031 


•151400 




IMC 


P,P 


0003P 


' 1 P64O0 




SUB 


1, 1 


00033 


•177777 




DVD 




00034 


' 131000 




MO V 


1,? 


0035 


•0O600PS 




,JSK 


B MPYU 


0003 6 


•050745 




STA 


P, SAVP 


00037 


' 1 1 1000 




MOV 


0,P 


00040 


' 0PO430 




LDA 


0, ATMCF 


00O41 


•0P44P6 




LDA 


1, ATCFCT 


00P/|P 


•(710 £,(710 1 J 




JSH 


©.POLY 



INT AKCTAM(OF TVO ARGUMENTS) 



J SAVE. RETURN ADDRESS 
% SAVE ACPI, 1 

; GET SIGM 

JIM ACP 

J SAVE IT 

J SAVE SIGM OF AC 1 

J AS SUFFLIMEMT FLAG 

J SET COMPLEMENT FLAG 
J IF AC0>AC1 

;and swap arguments 



JAC0<=AC1 

J CANNOT BE EOUAL 
JAC0/ACP < 1 

J GE'I AC0/ACP 

J GET X**P 



i PO I M T A R T A M CO F F F . S 
JARCTAN COEFF COUNT 
DEVALUATE POLYNOMIAL 
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PROGRAM LISTING: 



ARCTANGENT (cont'd) 



(Single Precision) 



00043*030740 

00044 '00 6002$ 

000 4 5 '01 4 740 
00046*000404 
00047 '02441 7 
■7100 50 ' 10 6400 
00051 • 121000 
00052' 101 PP0 
000 53 '014733 
00054*000404 
00055 * 105000 
000 5 6 '0204 10 
000 5 7 » i P2400 
00060 '024724 
000 6i ' l 25220 
000 62' 1 1 1200 
000 63 '02071 6 
000 64 '02471 6 
00065 '00271 3 

00066*144417 PI. 2: 
00067*00000 6 ATCFCT 
00070*000071 ' ATNCF: 
00071 '000403 
00072' 1 7550?. 
00073 '005275 
00074M 67274 
00075*014563 
00076' 1 52527 
00077 '077777 



LDA 

r> ( • '/ 
l.'O ~ 

JMP 

LDA 

SUR 

MOV 

MO VZ H 
n c*2 

JMP 

MO V 

LDA 

SUB 

LDA 

M0V2K 

MO Vk 

LDA 

LDA 

JMP 

14^417 

6 

. + 1 

402 

-2276 

5275 

- 10504 

145 63 

-25251 

77777 



2, SAV2 
9 MPYU 
COM PL 
.+4 

1,PI.2 
0, 1 

1,0 
0,0 
3 UP PL 

.+4 

0, 1 

0,PI. ? 
1,0 
1,3IGN 

1, 1 
0, 2 

0, SAV^ 
1,SAV1 
PKTURM 



.?GET X*P 

'COMPLEMENT ANGLE? 
J DO MOT COMPLEMENT 
', PI/2 
; Pi. P-x+-p 



i SUPPLIMENT? 

;no 



; PI -ANGLE 

i GET SIGN 
5 BIT IN 



SPI/2 1.44417 OCTAL 



J. 78 633 7 627 -2 
J-. 3701? 9999* -1 
J .83871 1896P -l 
; -. 1348 7 19 133 
; . 1988 1 48243 4 
;-, 33326 51491 7 
J. 99999 93478 2 



END 



1EMD OF ARC TAN 



BCD to BINARY 
(Single Precision) 

PURPOSE: 



TITLE: 



INPUT: 



OUTPUT: 



This routine converts a single precision number in BCD to 
binary. 

The title is . BCDB. 

A BCD integer in AC1 (maximum value 9999 decimal). 

Binary equivalent of BCD integer is returned in AC1. 



CALLING SEQUENCE AND ENTRY POINT: 

Indirect at . BCDB with normal return to the instruction 
following the call. 

ERROR CONDITIONS: 

If a digit greater than binary 1001 is encountered in the 
input, Carry will be set, AC1 will be unchanged, and ACO 
will contain the bad digit. Otherwise, Carry will be zero 
on return. 

CARRY AND REGISTERS: 

ACO, AC1, AC3, and Carry are destroyed; AC2 is unchanged. 

LENGTH AND TIME: 

This routine consists of 53 (octal) words and is normally 

relocatable. 

Execution time is 1.034 m.s. 
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BCD to BINARY 
(Double Precision) 



PURPOSE: 



TITLE: 



INPUT: 



OUTPUT: 



This routine converts a double precision number in BCD to 
binary. 



The title is . DBCB. 



A double precision integer is passed in ACO, AC1 (high order, 
low order) of maximum value of 99999999 decimal. 



The binary equivalent of the input is returned in ACO, AC1 
(high order, low order). 

CALLING SEQUENCE AND ENTRY POINT: 

Indirect at . DBCB with normal return to the instruction 

ERROR CONDITIONS: 

If a digit greater than 9 is encountered in the input, 
Carry will be set and ACO will contain the bad digit. 
Otherwise, Carry will be zero. 



CARRY AND REGISTERS: 

ACO, AC1, AC3, and Carry are destroyed; AC2 is unchanged. 

LENGTH AND TIME: 

This routine consists of 76 (octal) words and is normally- 
relocatable. 
Execution time is 2. 174 ms. 



BINARY TO BCD 



(Single Precision) 
PURPOSE: 



TITLE: 



INPUT: 



OUTPUT: 



This routine converts a binary number to its BCD equivalent. 



The title is . BBCD. 



An unsigned binary number in AC1. 



The BCD equivalent in AC1. 



CALLING SEQUENCE AND ENTRY POINT: 

Indirect to . BBCD with normal return to the instruction 
following the call. 

ERROR CONDITIONS: 

If a number greater than 9999 is input for conversion, no 
conversion will take place and Carry will be set. Otherwise, 
Carry will be zero. 

CARRY AND REGISTERS: 

AC1, AC3 and Carry are destroyed; ACO, and AC2 are unchanged. 

LENGTH AND TIME: 

This routine is 41 (octal) words and is normally relocatable. 
Execution time is 273. 8 -!- N * 14. 1 us where N is the 
sum of the digits of the result. 
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BINARY to BCD 



(Double Precision) 



PURPOSE: 



TITLE: 



INPUT: 



OUTPUT: 



This routine converts a double precision binary number to a 
BCD number. 



The title is . DBBC. 



A positive, double precision binary number in ACO, AC1 
(high order, low order). 



The BCD equivalent is in ACO, AC1 (high order, low order). 



CALLING SEQUENCE AND ENTRY POINT: 

Indirect to . DBBC with normal return to the instruction 
following the call. 

ERROR CONDITIONS: 

If AC 0, AC1 contains a number greater than 99999999, no 
conversion will take place and Carry will be set. Otherwise, 
Carry will be reset. 

CARRY AND REGISTERS: 

AC(\ AC1, AC3, and Carry are destroyed; AC2 is unchanged. 

LENGTH AND TIME: 

This routine consists of 57 (octal) words and is normally 
relocatable. 
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BINARY to DECIMAL 



(Single Precision) 
PURPOSE: 



TITLE: 



INPUT: 



OUTPUT: 



This routine converts a single precision two's complement 
number to an ASCII character string. 



The title is .BIND. 



A single precision, two's complement integer is passed in 
AC1. 



An ASCII character string terminated by a null word. 
Characters are passed right adjusted in ACO to the 
routine whose address must be in ZREL location .PTCH. 
The string is of the form: 

+DDDDD(NULL) 

-DDDDD(NULL) 



or 



CALLING SEQUENCE AND ENTRY POINT: 

Indirect at . BIND with normal return to the instruction 
following the call. 

CARRY AND REGISTERS: 

AC1, AC3 and Carry are destroyed; ACO and AC2 remain 
unchanged. 

LENGTH AND TIME: 

This routine consists of 51 (octal) words and is normally 

relocatable. 

Execution time is (378. 3 + N * 14. 1) ^s where 

N is the sum of the digits of the result. 
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BINARY to DECIMAL 



(Double Precision) 



PURPOSE: 



TITLE: 



INPUT: 



OUTPUT: 



This routine converts a double precision two's complement 
number to an ASCII decimal character string. 



The title is . DBD. 



A double precision, two's complement integer is passed 
in AC1, AC2 (high order, low order). 



ASCII character string of the form: 

or +DDDDDDDDDD(NULL) 

-DDDDDDDDDD(NULL) 
is outputted. Characters are passed right adjusted, 
bit 8 = 0, in ACo to a user routine whose address must 
be stored in ZREL location = PTCH. 

CALLING SEQUENCE AND ENTRY POINT: 

Indirect to . DBD with normal return to the instruction 
following the call. 

CARRY AND REGISTERS: 

AC1, AC2, AC3, and Carry are destroyed; ACO remains 
unchanged. 

LENGTH AND TIME: 

This routine consists of 112 (octal) words and is normally 

relocatable. 

Execution time is 1. 061 + N * . 047 ms where N is 

the sum of the digits of the result. 
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BINARY to GRAY CODE 



PURPOSE: 



TITLE: 



INPUT: 



OUTPUT: 



This routine computes the Gray Code equivalent of a 16-bit 
binary word. 



The title is . BGRY. 



A binary word in ACO. 



Gray Code equivalent in ACO. 



CALLING SEQUENCE AND ENTRY POINT: 

Indirect to . BGRY with normal return to the instruction 
following the call. 

CARRY AND REGISTERS: 

ACO, AC3, and Carry are destroyed; AC1, AC2 are unchanged. 

LENGTH AND TIME : 

This routine consists of 13 (octal) words and is normally 

relocatable. 

Execution time is 50. 3 u s. 
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BINARY to OCTAL 



(Single Precision) 



PURPOSE: 



TITLE: 



INPUT: 



uUiPUi: 



This routine converts a 16 -bit binary word to an octal 
ASCII character string. 



The title is . BINO. 



A 16-bit binary number is passed in AC1. 



An ASCII character string terminated by a null character. 
Characters are passed right adjusted in ACO to the user 
routine whose address must be stored in ZREL location .PTCH. 
The string is of the form: 

OOOOOO(NULL) 
where "O" represents octal digits. 

CALLING SEQUENCE AND ENTRY POINT: 

Indirect to . BINO with normal return to the instruction 
following the call. 

CARRY AND REGISTERS: 

ACO, AC1, AC3, and Carry are destroyed; AC2 remains 
unchanged. 

LENGTH AND TIME: 

This routine consists of 27 (octal) words and is normally 

relocatable. 

Execution time is 367. 6 FN* 20. 0/iswhere N is 

the sum of the digits of the result (the sum expressed in 

decimal). 
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DECIMAL TO BINARY 
(Single Precision) 

PURPOSE: 



TITLE: 



This routine converts ASCII characters to a single precision 
binary number. 

The title is . DBIN. 



INPUT: 



Input characters will be requested by calling a user "get a 
character" routine whose address must be stored in ZREL loca- 
tion . GTCH. This user routine must be provided. ASCII 
characters should be returned, right adjusted in ACOwith 
bit 8=0. This routine need not save any registers or Carry. 
Input should be in the form: 

S D D . . . D D (break) 
where "S" represents the sign ("-" or optionally "+"), D 
represents an ASCII decimal digit, and "break" is any ASCII 
character other than a digit. 

OUTPUT: 

Upon exit,AC0 will contain the ASCII break character and 
AC1 will contain the single precision, two's complement 
binary equivalent of the input. 

CALLING SEQUENCE AND ENTRY POINTS: 

Indirect to .DBIN with normal return to the instruction 
following the call. 

If it is desired to output a signal character, the calling sequence is 
indirect to . DBNI. An ASCII "S" followed by a null character 
will be transmitted via ACO to a "put a character" routine 
whose address must be in ZREL location .PTCH. 



ERROR CONDITIONS: 

Caution: The absolute value of the result is N MOD 2**15. 
For example: +96741 converts to +31205. 
-2**15 converts to 0. 

CARRY AND REGISTERS: 

ACO, AC1, and Carry are destroyed; AC2 and AC3 are unchanged. 

LENGTH AND TIME: 

This routine consists of 65 (octal) words and is normally relocatable. 
Execution time is approximately 110 f I * 82.2 us where I 
is the number of digits in the input. 
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DECIMAL to BINARY 



PURPOSE: 



TITLE: 



(Double Precision) 



This routine converts ASCII characters to a double precision 
binary number. 



The title is . DDB 



INPUT: 



Input characters will be requested by calling a user "get a 
character" routine whose address must be stored in ZREL location 
. GTCH. This user routine must be provided. ASCII 
characters should be returned, right adjusted in ACO with 
bit 8 = 0. This routine need not save any registers or Carry. 
Input should be in the form: 

"--''-'• • • — -— \"- — "/ 
where "S" represents the sign ("-" or optionally "+"), D rep- 
resents an ASCII decimal digit, and "break" is any ASCII 
character other than a digit. 



OUTPUT: 

Upon exit, ACO will contain the ASCII break character. AC1 
and AC2 contain the double precision two's complement equiva- 
lent of the input. 

CALLING SEQUENCE AND ENTRY POINTS: 

Indirect to . DDB with normal return to the instruction following 
the call. 

To output a signal character, the call is indirect to . DDBI . 
This will cause an ASCII "D" followed by a null to be sent via 
ACO to a "put a character" routine whose address must be 
stored in ZREL location .PTCH . This routine must accept 
ASCII characters in the same format as . GTCH . 

ERROR CONDITIONS: 

Caution: The absolute value of the result is N MOD 2**31 . 
(see .DBIN). 

CARRY AND REGISTERS: 

All accumulators and Carry are destroyed. 

LENGTH AND TIME: 

This routine consists of 77 (octal) words and is normally 
relocatable. Execution time is approximately 
69. 90 + I * 43. 35 us on the Nova 1200, where I 
is the number of digits input. 
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Signed DIVIDE 



(Single Precision) 
PURPOSE: 



TITLE: 



INPUT: 



OUTPUT: 



Divides two fixed point , two's complement numbers. 



The title is . DIV. 



Dividend in AGO (high order) and AC1 (low order). 
Divisor in AC2. 



Quotient in AC1. Remainder in ACQ (same sign as dividend). 



CALLING SEQUENCE AND ENTRY POINT: 

Indirect to . DIV with normal return to the instruction 
following the call. 

ERROR CONDITIONS: 

If the magnitude of the quotient exceeds 2**15-1, Carry is 
set and the dividend remains unchanged. Otherwise, Carry 
will be zero. 

CARRY AND REGISTERS: 

ACO, AC1, AC3, and Carry are destroyed; AC2 remains 
unchanged. 

LENGTH AND TIME: 

This routine consists of 45 (octal) words and is normally 

relocatable. 

Total average execution time is 605 ps. 
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Signed DIVIDE 
(Double Precision) 



PURPOSE: 



TITLE: 



liNnji: 



This routine calculates the quotient of two signed double 
precision numbers. 



The title is DDIV. 



The double precision divisor must be in ACO, AC1 (high order, 
low order). The quadruple precision dividend should be stored 
in four consecutive words, highest order to lowest order. AC2 
must contain the address of the highest order word of the dividend. 



OUTPUT: 



THP HrmWIp nrppicinn nnnHpnt io -rvM-n-*-T-»i=>H i-n AP/1 



a m /u-; ~U 



order, low order). Its sign is determined by the algebraic 
rules for signed division. The double precision remainder 
is stored in two consecutive memory words with the high order 
word first. AC2 will contain the address of the higher order 
remainder word. The sign of the remainder is the same as 
the sign of the dividend. 



CALLING SEQUENCE AND ENTRY POINT: 

Indirect to . DDIV with normal return to the instruction 
following the call. 

ERROR CONDITIONS: 

If the magnitude of the quotient would exceed 2**31-1 or| dividend 
>!divisorl, an error condition exists, Carry is set, and return 
is made with unpredictable results. 

CARRY AND REGISTERS: 

All accumulators and Carry are destroyed. 

LENGTH AND TIME: 

This routine occupies 140 (octal) locations. Execution 
time is approximately 2. 98 milliseconds. 

REFERENCE: 

"How to Use the Nova Computers ", section 2. 2. 
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Unsigned DIVIDE 



(Single Precision) 
PURPOSE: 



To calculate the quotient of two unsigned integers. 



TITLE: 



ENTRIES: 



INPUT: 



DVD 



DVD, .DIVI, .DIVU 



OUTPUT: 



If DVD or . DIVU entry, the dividend is input in ACO (high order 
portion) and AC1 (low order portion). If . DIVI entry, the 
dividend is input in AC1 alone. In either case, the divisor is 
input in AC2. 



The remainder is output in ACO, the quotient in ACL 



CALLING SEQUENCES: 

The calling sequences consist of indirect calls through 
page zero entries . DIVI or . DIVU with return to the 
next sequential location following the call. DVD is 
equivalent to JSR @ . DIVU. 

ERROR CONDITIONS: 

Any inputs which would yield a quotient larger than 2 - 1 
(i. e. , ACO > AC 2) causes Carry to be set and return to 
be made to the caller. Otherwise Carry is reset. 

CARRY AND REGISTERS: 

AC2 is unchanged; ACO, AC1 and AC3 are destroyed 
under normal operation. AC1 is also left unchanged 
under error conditions, and Carry is always set or 
reset as discussed above. 

LENGTH AND TIME: 

Two ZREL locations and 21 octal NREL locations are 
required. Average execution times for . DIVI are 103 m-s 
on the Supernova and 545 nson the Nova. Average exe- 
cution times for DVD are 102 [ison the Supernova and 
539 M-s on the Nova. 
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Unsigned DIVIDE (cont'd) 

(Single Precision) 

ALGORITHM: 

. DIVI creates an all-zero high order dividend upon entry 
and then enters the DVD common logic. 

DVD initially compares the divisor with the high order portion 
of the dividend. If that is greater than or equal to the divisor, 
the result would exceed 2 lS -i and could not be represented in 
16 bits usine r conventional twn's mmnlomont nntofinr t„ +ui„ 

u *■- -■• — — w.**.f .i.w**iv*lll. UVWAL1UU* 111 LlllO 

case Carry is set and return is made. 

Otherwise, the less significant (LS) half of the dividend is 
shifted left, and 16 iterations of the following logic are 

L - -*«. ^ ^w*. w*C uiuiv ox^iixxxucLiiL ^ivio; llctXl Ul Lilt 

dividend (ACO) is shifted to the left and is compared to the 
divisor (AC2). If the divisor is equal to or less than the 
MS dividend, the divisor is subtracted from the MS dividend 
and the LS dividend (AC1) is shifted left. Otherwise, no 
subtraction is performed and only the LS dividend left shift 
is performed. In both cases, Carry contains the latest 
quotient bit and is shifted in behind the MS portion of the 
dividend when the iteration is repeated. Upon completion, 
the 16-bit quotient is entirely assembled in AC1, and the 
final adjusted dividend in ACO is the remainder. The 
routine yields an exact answer. 



REFERENCE: 



"How to Use the Nnva rnmnm-prc" oo^nr.-n n n 
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Unsigned DIVIDE (cont'd) 
(Single Precision) 



00000-000001 • .DIVU: 
00001-000000* .DIVI: 



006000- 



. TITLE DVD 

.ENT DVD*. DIVI*. DIVU 

.ZREL 

DIVU 

DIVI 



.NREL 
DVD = JSR 



§.DIVU 



00000*102400 DIVI: SUB 0*0 

00001*054416 DIVU: STA 3*SAV3 

00002' 142432 SUBE# 2*0*S£C 

00003'000412 JMP DIVI 

00004*034414 LDA 3*M20 

00005*125120 M0V2L 1*1 

00006*101100 DVD0: MOVL 0*0 

00007*142412 SUB# 2*0*SZC 

00010' 142400 SUB 2*0 

00011*125100 MOVL 1*1 

00012*175404 INC 3*3*SZR 

00013'000773 JMP DVD0 

00014M76441 SUBO 3*3*SKP 

00015*176420 DIVI: SUBZ 3*3 

00016*002401 JMP @SA73 



INTEGER DIVIDE* CLEAR AC0 

SAVE AC 3 

TEST FOR OVERFLOW 

SET CARRY AND RETURN 

16 ITERATIONS 

SHIFT LOW DIVIDEND 

SHIFT HIGH DIVIDEND 

DOES DIVISOR GO IN? 

YES 

SHIFT LOW DIVIDEND 

CHECK COUNT 

NOT DONE 

DONE * CLEAR CARRY 

SET CARRY 

RETURN 



00017*000000 SAV3: 
00020' 177760 M20: 



; SAVE AC 3 

-20 ; - 16 DECIMAL 

.END ;END OF INTEGER DIVIDE 
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GRAY CODE to RTNARV 



PURPOSE: 

This routine computes the binary equivalent of a 16-bit 
Gray Code word. 

TITLE; 

The title is . GRYB. 

INPUT: 

Input is a Gray Code word in ACQ. 

OUTPUT: 

The binary equivalent is returned in AGO. 

CALLING SEQUENCE AND ENTRY POINTS: 

Indirect to . GRYB with normal return to the instruction 
following the call. 

CARRY AND REGISTERS: 

ACO, AC3 and Carry are destroyed; AC1, AC2 are unchanged. 

LENGTH AND TIME: 

This routine consists of 22 (octal) words and is normally 

relocatable. 

Execution time is 536.4 us. 
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Logical Exclusive OR 



PURPOSE: 

This routine computes the exclusive OR of two unsigned 
numbers. 

TITLE: 

The title is . XOR. 

INPUT: 

One 16 -bit quantity is passed in ACO, the second in AC1. 

OUTPUT: 

The exclusive OR of the two quantities is returned in ACO. 

CALLING SEQUENCE AND ENTRY POINT: 

Indirect to . XOR with normal return to the instruction 
following the call. 

CARRY AND REGISTERS: 

ACO, AC3, and Carry are destroyed; AC1 and AC2 are unchanged. 

LENGTH AND TIME: 

This routine is 7 words and is normally relocatable. 
Execution time is 34. us. 
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Logical Inclusive OR 



PURPOSE: 

This routine computes the logical inclusive OR of two 
unsigned numbers. 

TITLE: 

The title is . OR . 

INPUT: 

One 16 -bit quantity is passed in ACO, the second in AC1. 



OU ifui : 



The inclusive OR of the two quantities is returned in ACO. 

CALLING SEQUENCE AND ENTRY POINT: 

Indirect to . OR with normal return to the instruction 
following the call. 

CARRY AND REGISTERS: 

ACO is destroyed; AC1, AC2, AC3, and Carry are unchanged. 

LENGTH AND TIME: 

This routine consists of 5 words and is normally relocatable. 
Execution time is 25. 6 us. 
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Signed MULTIPLY 



(Single Precision) 
PURPOSE: 



TITLE: 



INPUT: 



OUTPUT: 



This routine multiplies two, fixed point, single precision, 
two's complement numbers. 



The title is . MPY. 



One fixed point, single precision operand is passed in AC1, 
the second in AC2. 



The double precision result is returned in ACO (high order) 
and AC1 (low order). 



CALLING SEQUENCE AND ENTRY POINT: 

Indirect to . MPY with normal return to the instruction 
following the call. 

CARRY AND REGISTERS: 

ACO, AC1, AC3, and Carry are destroyed, AC2 remains unchanged. 

LENGTH AND TIME: 

This routine consists of 16 (octal) words and is normally 

relocatable. 

For execution time in addition to unsigned multiply, 

56.4 us. 

For execution time in unsigned multiply, 340 us. 

Total average execution time is 396. 4 us. " , 
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Unsigned MULTIPLY 
(Single Precision) 



PURPOSE: 



TITLE: 



ENTRIES: 



INPUT: 



To calculate the product of two unsigned integers or the sum 
of that product and a third unsigned integer (MPY). 



MPY 



MPY, .MPYU, .MPYA 



OUTPUT: 



The multiplier and multiplicand are input in AC1 and AC2 (MPY). 
If a third integer is to be added to the product, the integer is input 
in ACO. 



The result is output with the high order portion in ACO and the 
low order portion in ACi. 



CALLING SEQUENCES: 

The calling sequences consist of indirect calls through page zero 
entries . MPYU and . MPYA, with return to the next sequential 
location following the call. MPY is equivalent to JSR @ . MPYA. 
. MPYU is used if no third integer is to be added to the product. 

ERROR CONDITIONS: 
None. 

CARRY AND REGISTERS: 

AC2 and Carry are restored; ACO, ACI, and AC3 are destroyed. 

LENGTH AND TIME: 

Two ZREL locations and 14 octal NREL locations are required. 
Average execution times for . MPYU are 87 us on the Supernova 
and 441 us on the Nova. Average execution times for MPY are 
86 us on the Supernova and 435 us on the Nova. 

ALGORITHM: 

. MPYU creates an all-zero addendum upon entry, and then enters 
the MPY common logic. The multiply -and -add function is performed 
by iteratively examining the least significant bit (LSB) of the multipli- 
cand and then shifting the multiplicand to the right. If the LSB was a 
one, then the multiplier is added to an accumulating partial product 
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Unsigned MULTIPLY (Continued) 

(Single Precision) 

ALGORITHM (cont'd) 

and that sum is shifted to the right. If the LSB was a zero, the 
accumulating partial product is simply shifted rightwards. 
Initially the partial product is zero in . MPYU, and is equal 
to the addendum in MPY. The process is carried out for 16^q 
iterations and yields the exact answer. 

REFERENCE: 

"How to Use the Nova Computers", section 2.2 . 
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Unsigned MULTIPLY (cont'd) 
(Single Precision) 



00000-000001 ' .MPYA; 
0000 1 -000000 • .MPYU : 



.TITLE MPY 

.ENT MPY , . MPYU,. MPYA 

.ZREL 

MPYA 
MPYU 



006000- 

ww www Iwcznow rir i u : 

00001' 0544 11 MPYA: 
00002*034411 
00003*125203 MPY1 : 
00004' 101201 
00005' 143220 
00006' 175404 
00007*000774 
00010' 125260 



0001 1 '002401 

00012'000000 SAV3: 
00013*177760 M20: 



.NREL 

MPY = JSR 

SUBC 0,0 

STA 3,SAV3 
LDA 3,M20 
MOVE 1,1, SNC 
MOVR 0,0, SKP 
ADDZR 2,0 
INC 3,3,SZR 
JMP MPY1 
MOVCR 1,1 

JMP @SAV3 


-20 



§ .MPYA 

I CLEAR ACQ, DON'T DISTURB 

; CARRY 

; SAVE AC 3 

; 16 TIMES THRU LOOP 

; CHECK NEXT MULTIPLIER BIT 

; 0, JUST SHIFT 

; 1, ADD MULTIPLICAND AND SHIFT 

; CHECK FOR 16TH TIME THRU 

i NO, CONTINUE 

; YES, SHIFT LAST LOW BIT 

; (NOTE IT WAS COMPLEMENTED BY 

,* FINAL INC) 

; RETURN 

; RETURN ADDRESS 

; -16 DECIMAL 



,END 



;END OF UNSIGNED MULTIPLY 
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Signed MULTIPLY 



(Double Precision) 



PURPOSE: 



TITLE: 



INPUT: 



This routine calculates the product of two signed double 
precision numbers. 



The title is DMPY. 



The double precision multiplier must be in ACO, AC1 (high 
order, low order). The double precision multiplicand should 
be stored in two consecutive words, higher order first. The 
location following the calling instruction (JSR) should contain 
the address of the high order word of the multiplicand. 



OUTPUT: 



The quadruple precision product is stored in four consecu- 
tive locations within the DMPY subroutine, highest order word 
first. AC2 contains the address of the highest order word. 
The sign of the product is determined by the algebraic rules 
for signed multiplication. 

CALLING SEQUENCE AND ENTRY POINT: 

Indirect to . DMPY with return to the second location after 
the call since the location after the call contains the address 
of the multiplicand. 

CARRY AND REGISTERS: 

All accumulators and Carry are destroyed. 

LENGTH AND TIME: 

This routine occupies 103 (octal) locations. Execution 
time is approximately 1. 62 milliseconds. 



REFERENCE: 



'How to use the Nova Computers", section 2. 2. 
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(Double Precision) 



PURPOSE: 



TITLE: 



INPUT: 



OUTPUT: 



This routine computes -D where D is a double precision 
two's complement integer. 



The title is . DNEG. 



A double precision, two's complement number in ACO, AC1 
(high order, low order). 



The negative of the input is returned in ACO (high order), 
AC1 (low order). 

CALLING SEQUENCE AND ENTRY POINT: 

Indirect to . DNEG with normal return to the instruction 
following the call. 

ERROR CONDITIONS: 

Caution: The negative of -2**31 cannot be represented and 
is returned unchanged. 

CARRY AND REGISTERS: 

ACO, AC1, AC3, and Carry are destroyed; AC2 remains unchanged. 

LENGTH AND TIME: 

This routine consists of 4 words and is normally relocatable. 
Execution time is 13. 8 us. 
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UCTAL to BINARY 



(Single Precision) 
PURPOSE: 



TITLE: 



INPUT: 



OUTPUT: 



This routine converts an ASCII octal character string to a 
binary number. 



The title is . OBIN 



Input characters will be requested by calling a user -written 
"get a character" routine whose address must be stored in 
ZREL location .GTCH . This user routine must be provided. 
Upon return from the call, this routine should return an 
ASCII character, right adjusted in ACO with bit 8=0. Input 
should be of the form: 

00. . . OO(break) 
where "O" represents octal digits. 



ACO contains the break character, and AC1 contains the 
binary number (MOD 200000 octal). 



CALLING SEQUENCE AND ENTRY POINTS: 
Indirect to . OBIN . 

It is desired to output a signal character, the calling 
sequence is indirect to . OBNI . 
An ASCII "0" followed by a null character will be 
transmitted via ACO to a user-written "put character" 
routine whose address must be stored in ZREL location 
. PTCH . In both cases, return is to the first word after 
the call. 

ERROR CONDITIONS: 

Caution: Result is N MOD 200000 (octal) e. g. , 576452* 
Converts to 176452. 

CARRY AND REGISTERS: 

ACO, AC1, AC3, and Carry are destroyed; AC2 is unchanged. 

LENGTH AND TIME: 

This routine consists of 42 (octal) words and is normally 

relocatable. 

Execution time for . OBIN is 63. + I * 70. 2 m s 

where I represents the number of digits in the input. 
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PARITY GENERATOR 



PURPOSE: 

This routine computes the even parity bit over a 16 -bit 
number and returns the bit in Carry. 

TITLE: 

The title is . PRTY. 

INPUT: 

A 16-bit number is passed in ACO. 

OUTPUT: 

The even parity bit over the contents of ACO will be returned 
in Carry. 

CALLING SEQUENCE AND ENTRY POINT: 

Indirect to . PRTY with return to the word following the call. 

CARRY AND REGISTERS: 

AC3 and Carry are destroyed; ACO, AC1, AC2 remain unchanged. 

LENGTH AND TIME: 

This routine consists of 16 (octal) words and is normally 

relocatable. 

Average execution time is 215. 4 us. 
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POLYNOMIAL EXPANSION 
(Single Precision) 

PURPOSE: 



TITLE : 



ENTRY: 



INPUT: 



OUTPUT: 



To calculate an integer Polynomial expansion series of the form P(x) 
P + P 1 x 1 + P 2 x 2 + . . . + P n x n 



POLYN 



.POLY 



The order of the polynomial is input as an integer in AC1. The argument 
x is a positive number input in AC1 with its binary point normally to the 
left of bit 0, the most significant bit. (There is no sign bit reserved 
since all inputs must be positive. ) Any shifting of the binary point 
is understood to shift the implicit point of the result in like fashion. 



16 bits of SIGNIFICANCE 



implicit binary point 

A pointer to the highest order coefficient is in ACO where all 
coefficients are two's complement numbers, usually in the following 
format: 




implicit binary point 



Any shifting of the binary point is understood to shift the implicit 
binary point of the result in like fashion. The coefficients are placed in 
a list with the following structure : 



n 

Vi 



increasing addresses 



The result, output in AC1, is in the same format as the coefficients, 
with shifting of the implicit binary point as required. 
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POLYNOMIAL EXPANSION (Continued ) 

(Single Precision) 

CALLING SEQUENCE: 

JSR indirect through page zero entry 

. POLY 

ERROR CONDITION: 
None. 

CARRY AND REGISTERS: 

AC2 is saved; ACO, AC1, AC3 and Carry are destroyed. 

T ninnmtT 1 X ti— t rrtTi r-r-< 

JL.I11NU 1 tl AiND 1 iivm : 

One ZREL location and 22 octal NREL locations. Execution time 
depends directly upon the order of the polynomial expansion. 

ALGORITHM: 

A cumulative partial product is formed by successively adding 
the next highest order coefficient not processed to a partial 
nroduct formed hv v and i-hp hiehe^t <~>rde">- ^<->^ffi<~iopt "/-.*■ n-™- 
cessed, and then by multiplying this partial sum by x. Only the 
most significant word of this product is retained. The process 
is performed iteratively, until all coefficients through P x 
have been processed. P Q is added to the final result. 

In equation form, 

P(x) + ((((P n x)+ P n -i)x). . . ) + P 
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POLYNOMIAL EXPANSION (Cont'd) 
(Single Precision) 



PROGRAM LISTING: 



; FIXED POINT POLYNOMIAL 

J INPUT: IN ACS 

; AS A ON F WORD INTEGER 

; COEFFICIENT POINTER IN AGO 

; COEFFICIENTS IN DESCENDING POWER 

; IN TWO'S COMPLEMENT FORM 

; ORDER OF POLYNOMIAL IN AC1 



5 OUTPUT: 



IN TWO 



COMPLEMENT 



FORM IN AC1 







.TITLE 


POLYN 






.ENT 


. POLY 






.EXTD 


MPYU 






.3REL 




0(?i000-0000f7ir* • 


. POLY : 


POLYN 

. N REL 




aWQflCA ' 05441 7 


POLYN: 


STA 


3, PL KIN 


(r,P}&.(x\ '040/11 7 




ST A 


r^, COEEF 


O0OOp 'f7i4441 7 




STA 


1, COUNT 


0O0O3 'OP 641 5 




LDA 


l,f>00EFF 


00004' 121 ISP 


LOOP: 


MO V^L 


1,0, soc 


00005* 1P4440 




NEGO 


1, 1 


000(716 '00 60(711 <5 




JSK 


o MPYU 


n MP! 7 ' 1 OSOf^p 




MOV 


0, 1,S:~0 


00010 ' 1 34400 




NEG 


1* 1 


'71';i0 1 1 ' a 1 40 7 




IS? 


COEEF 


0OO1P'OPP4O6 




LDA 


r-,^ COEFF 


OOP 13 ' 107"0w 




ADD 


0, 1 


PI 1 /i ' 1 ^40 5 




DSF: 


COUNT 


000 15*000 7 67 




JMP 


LOOP 


OOO 1 6*00^4^1 




JMP 


OPLETN 



SAVE RETURN 
■ S A VE CO E F E. PO I N TE K 

AND COUNT 
i HIGHEST COEFF 
: NEGATIVE? 

; YES, MAKE IT + VE, SET CARRY 
: ACO* 1 = A01*ACP 
: NEGATIVE SIGN? 
: YES, NEGATE RESULT 
: RAISE COEFF. POINTER 
:ADD NEXT COEFF 
;T0 RUNNING SUM 
(DONE? 

: GO TO NEXT COEEF 
: RETURN 



OOOl 7 '00000(71 PLRTN: o 
000OO '000000 COEEF: 
00qpi '000(7100 COUNT: 



.END 



J END OF POLYNOMIAL 
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RANDOM NUMBER GENERATOR 



PURPOSE: 



TITLE: 



INPUT: 



OUTPUT: 



This routine generates a (pseudo) random sequence of integers 
in the range to 2**16-1. 



The title is . RAND. 



The address of the previous random value(or initially a 
starting value) must be provided in the word after the call 
to . RAND. 



The new 16-bit random result will be returned in ACO and 
will also replace the previous value in memory. 



CALLING SEQUENCE AND ENTRY POINT: 

Indirect to . RAND followed by the address of the old 
value. Return will be to the instruction after the 
address parameter. 

ERROR CONDITIONS: 

Caution: If a K-bit number (1 < K « 16) is needed, use the most 
significant K bits (the least significant K bits are 
not as random). For example, to obtain random 
N MOD 2, use the sign bit of the result. 

CARRY AND REGISTERS: 

ACO, AC3, and Carry are destroyed; AC1 and AC2 are unchanged. 

LENGTH AND TIME: 

This routine consists of 36 (octal) words and is normally 

relocatable. 

Execution time is 244. 7 us. 
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SINE, COSINE 



(Single Precision) 
PURPOSE: 



TITLE: 



ENTRIES: 



INPUT: 



To calculate the sine or cosine of an angle expressed in 
radians. 



SINX 



SINX, .COSX 



The argument m is input in ACO, within the range 
-4«=n-< 4, in the following format: 



sign {integer j fraction j 



OUTPUT: 



bit bits 1 and 2 bits 3 through 15 
The sign bit is set to a 1 only if the argument is negative. 



The result is output in AC1, in the same format as the 
input. 



CALLING SEQUENCES: 

JSR indirect through either page zero entry . SINX 

or .COSX as appropriate. 

Return is to the next sequential location following the call. 

ERROR CONDITIONS: 
None. 

CARRY AND REGISTERS: 

ACO is saved; AC2, AC3, and Carry are destroyed. 

LENGTH AND TIME: 

Two ZREL and 65 octal NREL locations. 

Average execution time on the NOVA 1200 is .9 ms. 
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_SINE, COSINE (cont'd) 



(Single Precision) 
ALGORITHM: 

Upon entry to . COSX, m is subtracted from A/2, and 
the sine logic is performed on this difference. If the original 
argument input to . SINX is negative, the result will be equal 
to -SIN(m); negative arguments input to . COSX will be treated 
as positive arguments, since COS(-m)=COS (m). 

Upon entry to . SINX the product of m and 2/tt is calculated 
by means of a call to . MPYU. 

m._ aeon KJt. uic oiinu logic consists 01 tne tollowing 
expansion: 

SIN ( — ^ ) = m (P + m (Pj + m (P 2 + m P 3 ))) 

with the following polynomial coefficients: 

P Q = .15707 * 10 
P! = -.64589 

P 2 = +. 79434 * 10" 1 

P 3 = -.43330 * 10" 2 

REFERENCE: 

John F. Hart, "Computer Approximations", New York: 

John Wiley & Sons, Inc. , 1968; pages 116 and 117, INDEX 3300. 
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SINE, COSINE (cont'd) 



(Single Precision) 

PROGRAM LISTING: 



FIXED POIMT SIN/ CO: 







.TITLE 


SIMX 






. ENT 


. SI NX* • COSX 






. EXTD 


.POLY* MPYU 






.-'KFL 




00000- 


000001 ' . SI NX: 


SI NX 




0000 1 - 


000000' .COSX: 


COSX 

. N KFL 




00000 ' 


12640 1 COSX: 


SUB 


1, 1,SKF 


00001 « 


126520 SINX: 


SUB2L 


1, 1 


0000 P' 


044456 


STA 


1,SNFLG 


0(1003 ' 


40457 


ST A 


0,SAV0 


A ' 


054455 


STA 


3,RTURN 


00005 


1 1 1 1 20 


MO V0L 


0,2 


0000 6 


125005 


MO V 


1 , 1 , SN K 


00007 


101020 


MO V7, 


0,0 


00010 


1265 60 


SUPCL 


1, 1 


0001 1 


'044453 


STA 


1,SIGN 


000 1? 


'0 2.443 7 


LDA 


1,. 20VPI 


000 1 3 


• 00 6002S 


JSK 


MPYU 


000 1 4 


•125120 


MOVSL 


1, 1 


0001 5 


' 10 1 100 


MO VL 


0,0 


0001 6 


' 1 52 5 60 


SUBCL 


2,2 


0001 7 


' 1 251 20 


MOV"L 


1, 1 


000P0 


*101 100 


MO VL 


0,0 


00001 


' 151 100 


MOVL 


2,2 


00000 


•014436 


DS5 


SNFLG 


000 03 


•151400 


INC 


2,2 


00004 


•151222 


MO V0K 


2, 2,S":C 


00005 


• 100000 


COM 


0,0 


171 Pf, 


' 1 51 23 2 


MOV£R# 


0,2,SF.C 


i71 p 7 


•01043 5 


I S" 


SIGN 


00030 


' 1 1 1000 


MO V 


0,2 


H0031 


'05043? 


STA 


2, SAV2 


0003P 


' 10 5000 


MO V 


0, 1 


00033 


•00 6002S 


JSK 


§ M PY U 



', SET SIN/ COS FLAG 



', SAVE INPUT 

J SAVE RETURN 

J GET SIGN 

;MAKE IT +VE IF COS ENTRY 



5 SAVE SIGN 

5 GET 2/ PI 

; AC0, 1 = A0 1*AC2 

;PUSH OUT INTEGRA PART 



', ac2 is integer part 

;sine? 

', no, incremfm1 tntege 

5 odd integer part? 

J YES, COMPLIMEr.N'i FRACTION 
iODD INTEGER/2? 
>'YES CHANGE SIGN 
.: ER AC 'II ON IN AC 2 
', SAVE X 

JX**P 
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SIN/COS PROGRAM LISTING (cont'd): 



SINE, COSINE (cont'd) 
(Single Precision) 



00034' 1 1 100CT 
00035*02041 6 
00036*024414 
000 3 7' 00 600 IS 
00040 '030423 
00041 '00 600 2S 
00042* 1 51220 
000 43*143 220 
0^0 44'f-124 420 
00045' 125220 
00046* 105200 
■" 47 * 20 4 1 3 
000 50 '00 241 1 



MOV 


0,2 


LDA 


0,SMCOF 


LDA 


1,CFCNT 


JSR 


©.POLY 


LDA 


2, SAV2 




"rf ".'J" I U 


MOVER 


2,2 


A nrvz i 




LDA 


1>SIGN 


MO VS R 


1, 1 


MO VK 


0, 1 


LDA 


0>SAV0 


JMP 


6 K TURN 



5ACO POTNTS TO COEFF.S 
iOKDFR OF POLYNOMIAL 
> E VAL U AT E PO L YN M T A! 



^ GET SIGN IN CARRY 

J P 1 ISH I N TO HE S UL T 

J RESTORE INPUT ARGUMENT 

'> RETURN TO CALL+1 



00051 '121377 
00052*000003 
00053 '000054 ' 
00054' 1 77562 
00055*005053 
00056*126524 
00057 '044420 



• 20VPI 
CFCNT: 
SNCOF: 



121377 

3 

. + 1 

-21 6 
50 53 
-51254 
44420 



.; 2/PI C. 50574 6 OCTAL, TIEAKED) 



S -. 0043331 
', .079 4343 
; -. 645R92S 
i 1. 570791 CADD 



1 LATER) 



00060 '000000 

000 61 '000000 
000 62 '000000 

000 63 '000000 
000 64*000000 



SNFLG: 

RTURN: 

SAV0: 

SAV2: 

SIGN: 












END 



5 END OF SINX,COSX 
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SORT 



PURPOSE: 



This routine sorts a table of pairs of words so that the 
first words of the pairs are in ascending order. 



TITLE : 



The title is RSORT. 



INPUT: 



A table containing pairs of words; the first word of each 
pair, the key word, must be an unsigned integer. AC1 must 
contain the starting address of the table; AC2 must contain 
the ending address of the table. 



OUTPUT: 

A table of pairs of words, sorted so that the key words are 
in ascending order. The key word and its accompanying 
data word are unchanged by the sort. 

CALLING SEQUENCE AND ENTRY POINT: 

Direct to RSORT with normal return to the instruction 

following the call. 

CARRY AND REGISTERS: 

ACO, AC3, and Carry are destroyed; AC1 and AC2 are 
unchanged. 

LENGTH AND TIME: 

This routine consists of 127 (octal) words and is normally 
relocatable. Execution time is approximately . 3 seconds 
for a table containing 1000 pairs of words. 
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SORT (cont'd) 



PROGRAM LISTING: 



^Pl RSORT 



;ROuTl:\iE 10' SHUFFLE A i A PL E OF PAIRS OF WJFOJS 
; SO THE KIR ST t'.ORDS OF t'HE PA Ti-SC KEY i-ORDS) 

;ake tm ascending ordej<. 

;'KEY rtOKD-DAiA .-.ORD PAIRT^G LEFT jaiCHA ,\!GEC . 



J (LOAD I st AC1 biilVllNlG ADDRESS JF fARLE 

;AMl") ADDRESS OF fHE LASi wORD OF iARLE IN AC?) 

? ,JSF RSORT 

5 R ETuRN 

J T YP I CA L I T M E : .3 S EC F'J F 1 1 r<l ^ p a t H S 







• T T i L E 


F PORT 






. E'vT 


h S /-■ i 






• \REL 




■IPn^.T) 


"i?146"i F-sorr: 


L!)A 


'1 , . S i A A 


"ir^oil 


•*3 4T4 6P 


STA 


n, SiAK 


n^iq p 


"IP? 6^'l 


sup':R 


•i , .1 


nil^i 3 


' H'V,bc so-. I'l : 


T ,s - 


:.-.: i A * 


'innn/i 


">!5 64o5 


b i A 


3.»'~<:v 1'AK 


iiinD ^ 


' '"1 1 ^ 4 S 4 


TSP 


Si A* 


qffli'Til 6 


' "'t 5P^b3 


STA 


P,'^b lAK 


•1'T?I'"7|7 


' t 3 5P '1 P 


MOv 


1 , 3 


T:'"1"! 1 1 


'^P445^ 


LDA 


1 ,CP 


1 1 "' 1 1 


1 6 7 'i^T 


ADD 


3, 1 


q a q i p 


1 3P1 I P 


ADCL* 


1 > ?, SE-C 


ir*q i 3 


0^44-3 


OMP 


SORT 6 


HTT] 4 ' 


154447 


STA 


3* TE^P 


Pl'3<3 1 5 


'3P54P1 SO TIP: 


LDA 


1 ,*, 3 


fMi*l 6 


1 P341 b 


A M D ■? 


1 ,1, SMS 


"11^1 7 ' 


^0.^4 pi 


.J MP 


SOEI3 


fJl^^O^ ' 


3 04440 


LDA 


1»CP 


mci?i • 


1 3P4PP 


SUP 


1 y P 


f50^PP' 


TPSP^O 


LDA 


l.fl.P 


"! ^ !l 9 .V 


1 P341 4 


AMD* 


1 >P> SE-R 


Pl^IflPV 


P^'l41 5 


J'/iP 


SORT 4 


^ ,^| 91 o 5 • 


^4^143 7 


STA 


P * P I 1 


^PP6* 


PP1 4*1 P 


LDA 


^ * 1 # 3 


PTTP7 ' 


iPbrni 


LDA 


1 , , P 



5 SET JP SiACK POIMiER 

5 REi' ny j -i 

; SA vE SEYUR.'M 

J SA v E HIGH P I ;M T Eh 
JLJi'i POT.MiER' 



;SCAM COMPLETE? 

; YES 

J SAVE LOv. POIiMTER 

;lovv key bit sej? 

J.siO 

J DROP HIGH POINl'ER 

;HIGH KEY PIT SET? 

;YES,.\iO i.-vAP 

; SA vE COMPARISON Fill 

5 SWAP KEY WORD- DATA WORD PAIRS 
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SORT (cont'd) 



PROGRAM LISTING (cont'd): 



110 3 !7 I '^41 111 




STA 


1,1,2 




10031 '045410 




STA 


1 , , 3 




10^ 32' IP 1 411 




LDA 


0,1,3 




111 33 '025111 




LDA 


1 , 1 ,2 




11134*^41 111 




SfA 


1,1,2 




11135'145401 




STA 


1,1,3 




"ill 3 6* "1 PI 40 6 




LDA 


0, PIT 


J RECOVER COMPARISON rill 


111 37 '124423 


S0RF3: 


LDA 


1,C2 




1114 r, i • i 37 7n,7j 




ADD 


1 , 3 


; PUMP LJ.'- PUImTEK 


11041 ' 1 5641 4 


S0.RT4: 


SUB* 


2, 3, S£P 


5 SCAN OVER? 


■5)^4?'^^753 




J MP 


S0RT2 


5 NO, CONTINUE 


00143' 101 222 




M0V2R 


1,1, s^c 


jall bit positions covered? 


11144'11141 6 




J MP 


SORT 5 


JYES 


11145'1P441 6 




LDA 


1, TEMP 


; NO, NEXT BIT 


11146*^14735 




JSK 


S0RT1 




11147 '03P41P 




LDA 


2,1? STAK 




01151*114733 




JSR 


S0RT1 




0^51 ' HI 121 




MOV?L 


1,1, SKP 




(1W52' 101 100 


S0RT5: 


MOVL 


1,0 




01153'a26406 


S0KT6: 


LDA 


1,13 STAK 




01054*11 4415 




DS? 


STAK 




11155' 13 640 4 




LDA 


3,?STAK 




11156'H 4413 




DS? 


STAK 




11157*111 411 




JMP 


0,3 




111 61' 1111 65' 


' .STAK: 


STAK1 






111 61 '111111 


STAK: 









010 62*000002 


C2: 


2 






10063*000000 


TEMP: 


1 






111 64' 111111 


BIT: 


9 






11114? 


STAK1: 


• RLK 


42 










• END JE'.MD OF RADIX EXCHANGE SORT 
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r T»-i „■; ™*, ^ j pnitjnr n/-\/-\'-i-' 

(Single Precision) ~~ 



PURPOSE: 

TITLE: 
ENTRY: 
INPUT: 
OUTPUT: 



This routine calculates the square root of an unsigned single 
precision integer. 



ISQRT 

.ISQR 

An unsigned fixed point number is input in ACO. 



The square root of the input., truncated to the nearest 
integer, is output as an integer in AC1. Accuracy is to 
within one binary digit. 

CALLING SEQUENCE AND ENTRY POINT: 

Indirect through page zero entry . ISQR with return to 
the next sequential instruction following the call. 

ERROR CONDITIONS: 

There are no error conditions. All input values are 
acceptable as unsigned integers. 

CARRY AND REGISTERS: 

ACO, AC2, AC3, and Carry are destroyed. 

LENGTH AND TIME: 

This routine consists of 1 ZREL and 24 octal NREL 
locations. Average execution time is 82. 4 us on both 
the Supernova and Nova 800, and 153 /us on the Nova 1200. 
Average execution time on the Supernova SC is 60. 5 /is. 

ALGORITHM: 

One of three different algorithms is usually selected when 
extracting a square root: the sum of odd integers method, 
Newton's iterative method of successive approximations, 
and the longhand partial quotient procedure. The partial 
quotient procedure was selected for this subroutine. 

The partial quotient method of square root extraction 
is the same as the pencil-and-paper procedure commonly 
used to extract square roots of decimal numbers. In 
this method the number n whose root is to be extracted, 
is first paired off in groupings of two digits each 
(n Q m, etc.) 
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Unsigned SQUARE ROOT (cont'd) 
(Single Precision) 



ALGORITHM (cont'd) : 



V 



26 01 

f 1 

n o "1 



The largest integer root less than or equal to the 

exact root of no is selected, qg, and becomes the most 

significant digit of the extracted root, qg is squared 

and is subtracted from the first pair, yielding remainder tq. 



% 



-.126 01 
V 25 P^nj 




r Q is multiplied by 100 and is added to n p forming a new 
partial dividend, d} (the first partial dividend, d equaled ng, 
the first pair of digits), dj is divided by 20 * q + Qp where 
q l is the largest possible integer such that q± * (20*q Q + q^ 
* dt. The process is repeated with each partial dividend until 
all digits in n have been processed. 



(20*q ^ 




The subroutine algorithm is similar to that described above, 
but is modified to process 16 -bit unsigned binary numbers 
instead of real decimal numbers. The subroutine considers 
n as 8 pairs of binary digits and repeats its procedure 
8 times. 
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Unsigned SQUARE ROOT (cont'd) 
(Single Precision) 



ALGORITHM (cont'd): 



This algorithm was selected since its use results in the most 
efficient coding and fastest execution time for the complete 
range of permissible input values for all Nova family 
machines. 

The sum of odd integers method given in "How to Use the 
Nova Computers" yields more efficient coding and faster 
execution times when input arguments are no greater than 

1000,^ Over the whfllp Tflricrp nf nprmiooiWo i,-,™, -t-,-, 

however, its average execution time is 412 jis on the 
Supernova. 

Newton's method requires either the use of hardware 

v^.^v, ^j. ciiv- ouuvvait uiviuc OUUJLUU Lillet), cUlU ItS 

average execution time over the entire permissible 
range of inputs is 114 us on the Supernova (with 
hardware divide). 

REFERENCE: 

Section 2.2 of "How to Use the Nova Computers " contains 
a discussion of the sum of odd integers method of square 
root extraction. 
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Unsigned SQUARE ROOT (cont* d) 



(Single Precision) 



PROGRAM LISTING: 



00 1 I SORT 



{SQUARE KOJ 1 UK AN UNSIGNED IN 1 f_ li fcl K 

{INIEGEk ARGUMEN'i IN AC0 

{OUTPUT: SQUARE ROOT TRUNCATED TO NEA.KES 1 

{INTEGER RETOrtNED IN AC1 

{REGISTERS AND CARRY DESTROYED 

5 CALLING SEQUENCE: 

.JSR a. I SORT 

REIU^N LOCATION 



00000-000002 * . i sqr: 



0000 0' 000000 
00 001 '000000 
00002" 0547 7 6 
00003*024420 
0000 4*0447 7 5 
000105' 126400 
0000 6' 152 400 



RTRN: 
CO UNI 
I 50 Rl 



• 2REL 

I SQRT 

.NREL 





ST A 

LDA 

ST A 

SUB 

SuB 



. 1 I TL E 

• EN 1 



I SORT 
• I SQR 



JRELUmM ADDRESS 

; COUNTING LOCATION 



3,RTrN 
1,C* 

1 , CO UNI 
1, 1 

2,2 



{SAVE RETURN ADDRESS 
{SET UP CO UNI * 

{ CLEAR PARTIAL QUOTIENT 
5 AND REMAINDER 



00007' 101 120 
00010' 151 100 



0001 1 ' 
00012 l 
0001 3 ! 
00014' 
0001 5 1 
00016' 



1 1 1 20 
151 100 
135120 
175140 
1 72432 
17 2420 



0001 7 ' 1251 00 
00020 '01 47 61 
00021 '0007 66 
00022*002756 



ISQl: M0V6L 0,0 J ADD (NEXT) HIGHEST 

MOVL 2,2 J TWO BITS OF ARG. 10 

MOV?L 0,0 J FOUR TIMES PAR II AL REMAINDER 

MOVL 2,2 i TO GET NEW PARTIAL DIVIDEND 

MUv/ZL 1,3 J NEW QUOTIENT IS TWICE THE OLD ONE 

MOVOL 3,3 J IF TWICE NEw QOOTIENT+1 (DIVISOR) 

SOB£# 3,2,SZC US EQUAL 10 OR LESS 1HAN 

SUBZ 3,2 {THE DIVIDEND, INCREMENT 

MOVL 1,1 {QUOTIENT AND SUBTRACT 

DSZ COUNT J DIVISOR FROM DIVIDEND 

JMP I SOI {10 GET NEW REMAINDER 

JMP f?RTRN {AND REPEAT 7 MORE TIMES 



0002 3 '0000 10 C* 



10 {EIGHT 

•END {END Of SQRT( INTEGER) 



0002 I SORT 



C8 


000023' 


1/19 


1/37 


COUNT 


000001 ' 


1/17 


1/20 


ISQl 


000007' 


1/24 


1/34 


ISQR'I 


000002* 


1/14 


1/1R 


RTRN 


000000* 


1/1 6 


1/18 


.ISQR 


000000- 


1/1 4 





1/33 



1/35 
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Unsigned SQUARE ROOT 
(Double Precision ) 



PURPOSE: 



TITLE: 



INPUT: 



This routine calculates the square root of an unsigned 
double precision integer, and expresses the result as a 
truncated single precision integer. 



DISQR 



An unsigned double precision fixed point number is input 
in ACO and AC1 (more significant half in Arm 



The square root of the input, truncated to the nearest integer, 
is output as an integer in AC2. Accuracy is to within one 

hina y\T Hi rrif- 

— —-—--,/ «^t,^v. 

CALLING SEQUENCE AND ENTRY POINT: 

Indirect through page zero entry . DISQ with return to the 
next sequential instruction following the call. 

ERROR CONDITIONS: 

There are no error conditions. All input values are 
acceptable as unsigned integers. 

CARRY AND REGISTERS: 

ACO, AC1, AC3, and Carry are destroyed. 

LENGTH AND TIME: 

This routine consists of 1 ZR.EL and 16 octal NREL 
locations. Average execution time on the Nova 1200 with 
software divide is 6 ms. 

ALGORITHM: 

Initially a trial square root of 177777 octal is assumed. An 
integer division of the input argument by this trial square root 
is then performed, and the arithmetic mean of this quotient and 
the trial root is calculated, yielding a new trial root. The 
process of dividing the new trial root by the input argument 
continues until the new root differs by less than two from the 
next most recent trial root. 
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Unsigned SQUARE ROOT (cont'd) 



(Double Precision) 



PROGRAM LISTING: 



oooi riser 



;SCUAPE FOOT OF AN UNSICNED TEC VCPD INTEGER 
;iiJPU T INTEGER Ii- ; AGO, 1 

; RESULT (TRUNCATED TO NEAREST ) INTEGER II-: ACS 
; ACCUMULATORS AN D CAPPY LOFT 

; CALLING SEQUENCE: 



(LOAD INTEGER IN ACO AND AC1) 

JSP g. FIFO 
PESULT LOCATION 



;pance of afgunent to (2**16-1)**; 
;no efpop meffagef 
;suppcfting pcutine fit- 



;< UNSIGNED SINGLE PPECISION 



iTEGEP DIVIDE) 



00000-000003' 

00000 '000000 
00001 'OOCOCO 
00002' 000000 

00003 '054777 
00004*040774 
00005 '044774 
00006* 152000 

OOOC-7'020771 
OOCIO'024771 
00C1 1 ' 177777 
00012' 133220 
00013 * 146654 
00014'000773 
00015*002765 







.TITLE 


DI SOP 






. EN" 


• DISC 






. EXTN 


DVD 




• ZFEL 






.DISC: 


DISOP 
.NFEL 






INTGO: 









INTGl't 









FT P.N : 









DI SOP: 


STA 


3 j PTPf'J 






STA 


CI NT GO 






STA 


I, INTG1 






ADC 


2.j c 




DISCI : 


LDA 


Qj INTGO 






L DA 


1, IMTG1 






DVD 








ADDZP 


1,2 






SUBOR# 


2, 1,SZP 






JNP 


DISP1 






J MP 


©FT FN 





; SAVE PETUPN ADDFESS 

J SAVE HI GH VCFD 

;and LOl." VCPD 

JTFIAL PGCT SET 1 77777 C OCTAL ) 



;lgad apgt t nent 
j divide it ey tpial foot 
;take the mean of quotient and 
; tpial foot to improve hoot 
;if nev foot differs by less 



', THAN TUG 



POM QUOTIENT, RETURN 



. EN D 



END OF SPUAPE FOOT OF INTEGE! 
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SUBTRACTION 



(Double Precision) 



PURPOSE: 



TITLE: 



INPUT: 



OUTPUT: 



This routine computes the difference of two double precision 
two's complement integers. 



The title is . DSUB. 



The minuend is passed in ACO, AC1, (high order, low order). 
The subtrahend must be in two consecutive memory words, 

UirrUnf n-~Ars~ t^\'\m,^A Utt 1 „,,,„ ». „„J„„ TU „ J JT_11 ! 

lllgHCl U1UC1 AU11UHCU JJJ JLUWCi U1UC1, 4.11C WU1U lUllUWlIlg 

the JSR should contain the address of the higher order 
word of the subtrahend. 



ine aouoie precision aitierence is returned in A(J0, AC1 
(high order, low order). 



CALLING SEQUENCE AND ENTRY POINT: 

Indirect to . DSUB followed by the address of higher order 
word of subtrahend. Return will be to the instruction 
following the address parameter. 

ERROR CONDITIONS: 

Caution: No check is made for overflow. 

CARRY AND REGISTERS: 

ACO, AC1, AC3, and Carry are destroyed; AC2 remains 
unchanged. 

LENGTH AND TIME: 

This routine is 15 (octal) words and is normally relocatable. 
Execution time is 54. 9 u s « 
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rjJrojJDJ^ 2 T ° REVBION 3 " F ™* RHLOCATABLE M ATO 



a^nt tc!utr S *" d6SCribed '" *" t0ll ° Wlng ""• T ^P h >"< "actions 



0^, i r\ 



« 8 « ^ Mgnea double precision division is now performed by the routine DMPY. 

Page 30 Signed double precision multiplication is now performed by the routine 
DDIV. The descriptions of DDIV and DMPY (page 19) replace that of 
. DPMD. 

Page 42 A description of RSORT, which sorts pairs of words in ascending order, 
has been added. 

Some character conversion routines now require a "get a character" 
and/or "put a character" routine whose address must be stored in ZREL 
locations . GTCH and . PTCH, respectively. This change is reflected in 
the following pages: 12, 13, 15, 16, and 17. 
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